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ly-^ , Abstract 

o, 

^O , We consider a particle that is subject to a constant force and scatters inelastically 

on a vibrating periodically corrugated floor. At small friction and small radius of the 
f^ ' circular scatterers the dynamics is dominated by resonances forming spiral struc- 

tures in phase space. These spiral modes lead to pronounced maxima and minima 
in the diffusion coefficient as a function of the vibration frequency, as is shown in 
C ' computer simulations. Our theoretical predictions may be verified experimentally 

K^ . by studying transport of granular matter on vibratory conveyors. 
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1 Introduction 



A ball bouncing inelastically on an oscillating plate under the action of a gravi- 
tational field appears to be a very simple dynamical system. Surprisingly, such 
a ball exhibits an extremely rich dynamics thus providing a prominent example 
for complexity in a seemingly trivial nonlinear dynamical system: Experiments 
indicated period-doubling bifurcations into chaotic motion [1,2,3,4] while the- 
oretical analyses classified different types of dynamical behavior in terms of 
phase locking, chattering, and strange attractors [5,6,7]. 

Oscillating surfaces, on the other hand, are often used in the field of granular 
material in order to drive a gas of granular particles into a nonequilibrium 
steady state [8,9,10]. A model that is somewhat intermediate between an in- 
teracting many-particle system and a single bouncing ball is the bouncing ball 
billiard (BBB) [11]. Here the surface is not fiat but periodically corrugated 
mimicking a fixed periodic grid of steel balls [10]. A single particle bouncing 
inelastically on this vibrating surface thus performs diffusive motion. Study- 
ing the BBB at large friction and large radius of the scatterers, the diffusion 
coefficient turned out to be a highly irregular function of the driving fre- 
quency [11]. As the dynamical reason for the largest maxima some integer 
frequency locking between bouncing ball and vibrating plate could be identi- 
fied, a phenomenon that is well-known for a ball bouncing on a flat vibrating 
plate [2,3,4,5,6,12,13]. Additionally, we detected an irregular structure on fine 
scales that is supposedly due to some subtle effects hke pruning of orbits under 
parameter variation [14]. 

In this Letter we show that yet there exists another microscopic mechanism 
creating non-monotonicities in the frequency-dependent diffusion coefficient. 
Our starting point is to study diffusion in the BBB at smaller friction and 
smaller scatterer radius than in Ref. [11]. In this situation we find a spiraling 
bouncing ball mode that changes under parameter variation. We show that this 
mode is again responsible for the emergence of local maxima and minima in the 
diffusion coefficient as a function of the vibration frequency. We conclude with 
a brief outlook towards consequences of our work for the standard bouncing 
ball as well as for transport on vibro-transporters. 



2 The bouncing ball billiard 



The model is depicted in Fig. 1: It consists of a floor oscillating with yf = 
—A sin(27r/t), where A = 0.01 and / are the amplitude respectively the 
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Fig. 1. Sketch of the bouncing bah bilhard: A point particle cohides inelastically 
with circular scatterers forming a periodic lattice on a line. Parallel to y there is 
an external field with constant acceleration g. The corrugated floor oscillates with 
amplitude A and frequency /. 

frequency of the vibration. ^ This floor is equipped with a periodic grid of 
circular scatterers of radius R whose centers are a distance of S* = 2 apart 
from each other. We now consider a point particle^ performing a free flight 
between two collisions in a gravitational field with acceleration g = 9800 || y. 
The particle's coordinates (a^^+i, 2/n+i) ^^^ velocities {Vx,n+i^ "^y^n+i) ^^ time 
tn+i immediately before the {n + l)th collision and its coordinates (a;^,?/^) 



and velocities (f^„, v^^) at time t^ 
related by the equations [11] 



immediately after the nth collision are 



Vn+l = yn+ V^,n{^n+l " tn) " g{tn+l " t„) 72 



■^x,n+l 



f. 



y,n+l 



'-'y,n " QV^n+l — ^n) • 



(1) 

(2) 
(3) 
(4) 



At a collision the velocities change according to 






(5) 
(6) 



where Vf is the velocity of the corrugated floor. The local coordinate system 
for the velocities v±_ and v\\ is given in Fig. 1. We assume that the scattering 
process is inelastic by introducing the two restitution coefficients a and (3 
perpendicular, respectively tangential to the surface at the scattering point. 

In Ref. [11] we studied the BBB for i? = 25, a = 0.5 and /5 = 0.99. Here we 
consider the case of i? = 15 and a = 0.7, which is closer to the experiments of 



^ We work with reduced units, where all quantities are dimensionless in terms of 
the unit amplitude ^o = ^rnm, respectively the unit frequency /o = IHz. 
^ Note that the radii of the moving particle and of a scatterer are additive, hence 
we only neglect any rotational energy. 




Fig. 2. The diffusion coefficient D{f) of the bouncing ball billiard Fig. 1 at friction 
a = 0.7 and scatterer radius i? = 15 as a function of the vibration frequency / 
of the corrugated floor. The graph consists of 124 data points with error bars at 
/ = 55, 70.8 and 82.1. The labels (a) to (e) refer to the respective phase space plots 
of Fig. 3. 

Urbach et al. [10]. In contrast to naive expectations, this simple variation of 
parameters profoundly changes the diffusive dynamics of the BBB, as we will 
show in the following. 



3 Deterministic diffusion, frequency locking, and spiral modes 



The diffusion coefficient D was computed from simulations according to the 
Einstein formula 



D 



lim 



< (X(t)-X(0))2 > 



2t 



(7) 



where the brackets denote an ensemble average over moving particles; for fur- 
ther numerical details cf. Ref. [11]. Fig. 2 shows that D exhibits local maxima 
and minima as a function of the driving frequency /. However, already at a 
first view this curve looks very different from the corresponding one in Ref. 
[11], where the diffusion coefficient was computed at different parameter val- 
ues. 
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Fig. 3. Projections of the phase space of the bouncing ball billiard Eqs. (1) to (6) 
onto positions and velocities along the y-axis at the collisions. The frequencies of 
the vibrating plate are / = (a) 54, (b) 58, (c) 60, (d) 66, (e) 72, (f) 78. The spiraling 
lines represent the analytical approximation Eqs. (9), (10). 



The major irregularities of Fig. 2 can be understood by means of the phase 
space projections displayed in Fig. 3. Here we have plotted the positions y and 
the velocities Vy of a moving particle immediately after its collisions with the 
plate. The onset of diffusion around / = 50 is characterized by the existence 
of many creeping orbits, where the particle performs long sequences of tiny 
little jumps along the surface hence moving parallel to y like the harmonically 
oscillating plate. This dynamics is reminiscent in Fig. 3 in form of incomplete 
ellipses around wj" = thus representing a simple harmonic oscillator mode. 

The first maximum of D{f) around / = 54 reflects a 1/1 frequency locking 
between particle and plate, as is shown in Fig. 3 (a) as the spot around Vy ^ 
100. In other words, the time of flight Tp of the particle is identical with 
the vibration period Tf of the floor. Additionally, the system may exhibit 
a dynamics where this behavior is interrupted by periods of creepy motion. 
Whether this happens depends on the initial conditions, hence the dynamics 
is non-ergodic. However, only for the peak around / = 55 we could detect 
such difficulties, whereas for / > 58 our numerics indicated ergodic motion. 
Frequency locking has been widely discussed for the simple bouncing ball 
[2,3,4,5,6,12,13] and has already been found to enhance diffusion in the BBB 
[11], where all the major peaks of the frequency-dependent diffusion coefficient 
could be identified in terms of integer frequency locking. 

Around / = 58 this resonance is completely destroyed resulting in a local min- 
imum of D{f). In Fig. 3 (b) this is represented by a "smearing out" of the 1/1 
frequency locking spot together with a dominance of creeps. Surprisingly, the 
next local maximum at / = 60 corresponds to a new structure in phase space: 
The granular particle locks into a virtual harmonic oscillator mode (VHO) sit- 
uated around Vy ~ 70, cf. Fig. 3 (c), which reflects the harmonic driving and 
again enhances diffusion. Around / = 62 this mode again becomes unstable 
leading to a local minimum. At / = 66 the VHO starts to "spiral out", see 
Fig. 3 (d), eventually resulting in a two-loop spiral. Fig. 3 (e). This causes 
a drastic increase of -D(/) on a global scale. Note that we could not detect 
any simple periodic motion on this spiral, hence this scenario may not be un- 
derstood as a simple period-doubling bifurcation. We furthermore emphasize 
that parallel to x the motion is always highly irregular, as is reminiscent in the 
existence of a diffusion coefficient. At / = 76 another loop develops, see Fig. 
3 (f) with / = 78 for a slightly advanced stage, before this structure settles 
into a three-loop spiral that explains the washed-out local maximum in D{f) 
around / = 80. A fourth loop emerges at / = 86 being completed around 
/ = 90, followed by the onset of a fifth one at / = 96. 

In order to understand the frequency dependence of this spiral more quanti- 
tatively we analyze the data presented in Fig. 3 using a simple argument of 
Warr et al. [6]: The time of flight Tp between two collisions of the particle with 



the plate is determined by Tp = 2Vy/g, hence one may calculate 

^:^^^^ (8) 

Tf 9 

and check whether this fraction yields some rational number indicating some 
frequency locking. Extracting Vy from Fig. 3 one can easily verify that the 
black spot in (a) indeed corresponds to /c ^ 1. Complete spiral loops at higher 
frequencies appear to be situated around fc + c A/, c G No, with fc — 60 and 
A/ ^ 10. Let Vy be the maximal values of the spirals in Fig. 3. We then find 
that k ~ 1.47 around / = 60, /c ^ 2.5 around / = 72, /c ~ 3.43 around / = 80, 
and fc ~ 4.5 around / = 90. The spiral mode thus locks into "smeared-out" 
half-integer resonances starting around fi + i A/, i e Nq, /« — 66. 

We remark that, for a ball bouncing on a flat plate, signs of such spiral modes 
were already observed experimentally [6] and described theoretically [3,4,5,7]. 
Particularly, Luck and Mehta [5] provided a simple analytical approximation 
for the coarse functional form of this mode. Their derivation can straightfor- 
wardly be repeated for Eqs. (1) to (6): Let us approximate the surface to be 
flat, and let us assume that there are no correlations between the collisions. 
By using Eqs. (2) and (5) we then arrive at 



y+ = -Asin(27r/ti) (9) 

77+1 = ag/2{ti - to) - A27r/(1 + a) cos(27r/ti) , (10) 

cp. Eq. (10) to Eq. (3.20) of Ref. [5], where to is the initial time at which the 
particle was launched from the plate and ti is the time of the first collision. For 
to we calculated the time at which the plate moves with maximum positive 
velocity, since here, for high enough frequency, a particle will be launched 
which previously is at rest. For ti we obtained the distributions of collision 
times tc from simulations yielding a range of < tc < tmaxif)- In order to 
ensure that ti — to > we then defined ti := tc + to. Consequently, Eqs. (9), 
(10) contain no fit parameter. Results from Eqs. (9), (10) are displayed in 
Fig. 3 (d) to (f). Naturally, this approximation does not work well for small 
collision times, respectively in the regime of creeping orbits. However, starting 
from the VHO it reproduces the whole structure of the spiral mode very well. 



4 Outlook and conclusions 



In this Letter we have studied the BBB for smaller friction and smaller scat- 
terer radius than in our previous work Ref. [11]. Computer simulations demon- 
strated that, again, the diffusion coefficient is a highly irregular function of 



the vibration frequency. However, for the present parameters only the first 
local maximuni could be identified in terms of a simple frequency locking. We 
showed that there exist local extrema due to what we called a spiral mode 
in the BBB. Whenever the vertical velocity of the particle, respectively the 
vibration frequency of the plate, is large enough such that the particle can lock 
into a multiple half-integer frequency of the plate, the spiral mode develops 
another loop. The coarse functional form of this spiral was extracted from the 
equations of motion by neglecting the radius of the scatterers and by assuming 
memory loss at the collisions. 

In Ref. [7] these spiral structures have been denoted as "chattering bands" 
of the simple bouncing ball, and it has been argued that they are rather 
stable against random perturbations. In the BBB the scatterers are defocusing 
inducing an additional deterministic chaotic component into the dynamics 
[14]. Interpreting the geometry of the BBB as a "perturbation" thus leads to 
the conclusion that spiral modes form states of the bouncing ball dynamics 
which are much more stable than simple frequency locking resonances. 

Interestingly, bouncing ball-type models have already been used since quite 
some time in order to understand transport on vibro-transporters, which are 
common carriers of agricultural material and related matter. For a corrugated 
version that very much reminds of the BBB and is widely used in industrial 
applications see, e.g., Ref. [15]. In order to transport particles on such devices 
one generates an average drift velocity by means of some symmetry break- 
ing. For tilted systems with a flat surface Hongler et al. [12,13] predicted an 
integer frequency locking leading to local extrema in the transport rate as a 
function of a control parameter, as was recently verified in experiments [16]. 
Grochowski et al. constructed a circular vibratory conveyor, where the trough 
was driven by asymmetric oscillations [17,18,19]. Both for a tracer particle in 
a layer of granular matter and for a single highly inelastic test particle they 
observed a non-monotonic current reversal of the transport velocity under vari- 
ation of control parameters. This ratchet-like effect was recently reproduced 
in simulations [20] again indicating some underlying frequency locking. 

Our work leads to the conclusion that in vibro-transporters not only the cur- 
rent, but also other transport property quantifiers such as the diffusion coeffi- 
cient may be irregular functions of control parameters. So far we have identified 
three basic mechanisms generating irregular parameter dependencies, which 
are (i) frequency locking [11], (ii) the spiral modes discussed here, and (iii) 
pruning effects [14] leading to fractal-like irregularities on finer and finer scales. 
Looking for irregular transport coefficients in experiments may thus pave the 
way for a very fine-tuned control of the dynamics of vibro-transporters. This 
may eventually provide a basis for practical applications such as, for example, 
sophisticated methods of particle separation [9] . 
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